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Abstract 



A special asymptotic solution of the Painleve-2 equation with small parameter is stu- 
died. This solution has a critical point t* corresponding to a bifurcation phenomenon. 
When t < t» the constructed solution varies slowly and when t > t* the solution 
oscillates very fast. We investigate the transitional layer in detail and obtain a smooth 
asymptotic solution, using a sequence of scaling and matching procedures. 



> 

: 1 Introduction 

m 
o 



In this work a special asymptotic solution for the equation Painleve-2 

e 2 u" + 2u 3 + tu = 1 (1) 



O 

o 

^ . is constructed as e — > 0. 

The behaviour of wanted solution differs in different intervals of the parameter t. The 
, -Ch ' qualitative behaviour of numerical solution is indicated in the figure (at e 2 = 0.1). Calcu- 

lations for another small values of e give pictures like this. Let us explain these numeric 
^ ', results using asymptotic theory for small e. 

In the area I the special solution is approximated by an asymptotic solution in which 
the leading term is a least root of a cubic equation 2u 3 + tu = 1. Corrections are algebraic 
functions of t. Such algebraic asymptotics becomes invalid near the point t*, where two 
of the roots of the cubic equation coalesce. In the left neighborhood of the point t* this 
asymptotics goes into an asymptotics of a special solution of the Painleve-1 equation with 
respect to new scaling variable r = (t — t*)^ -4 / 5 . This is the area II in the figure. This 
special solution of the Painleve-1 equation has poles at r = r*., k = 0, 1, . . . . In the 
neighborhoods of these poles one more scaling is done. New variable is 9 = (r — T/^e -1 / 5 
(area III). In this area the asymptotics is defined by separatrix solution of a nonlinear 
autonomous equation. This combined asymptotic structure becomes invalid as r — ► oo, 
because the poles of the solution Painleve-1 equations close to each other. As r — > oo a 
fast oscillating asymptotics is valid. It is the area IV in the figure. 

The qualitative behaviour of solutions of second-order ordinary differential equations 
with respect to an additional parameter was explained, for example, in the book p|. In jfl] 
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the various types of bifurcations for equilibrium positions of conservative second-order 
ordinary differential equations were described also. 

The equation ([!]) is nonautonomous, however one can separate the dependence on slow t 
and fast variables (r, 6) in the asymptotic solution of this equation. We may speak about 
the phase plane and phase trajectory of this equation with respect to the fast variable t/e. 
The equilibrium positions on the phase plane depend on t. When t < t* there are three 
equilibrium positions. At the critical point i* the bifurcation "saddle-center" occurs. It 
means that one of stable and unstable equilibrium positions coalesce. When t > i* only 
one equilibrium position exists. Such bifurcation leads to instability 0]. 

The bifurcations of slowly varying equilibrium positions of second order equation with 
the algebraic nonlinearity and slowly varying parameters were considered in || only in 
a preliminary fashion. When this text was written, new work [Q] became known. In that 
work, a change of an energy and a phase jump has been studied for a solution in a very 
narrow layer near a saddle-center bifurcation point in general case. However results of the 
work Q are inapplicable to the Painleve-2 equation, because this equation has degenerate 
behaviour with respect to equations considered in 0. 

Amongst another works, in which the asymptotic solutions of the nonlinear equations 
with varying coefficients were studied we should note the work |5|. In that work a chan- 
ging of an adiabatic invariant was studied in a problem when a solution passed through 
separatrix in the nondegenerate case. One more work where the passage through the 
separatrix in nondegenerate case was studied is Q. 

More complicated bifurcation is the pitchfork in the equation Painleve-2. Asymptotics 
with respect to a small parameter of the solutions for the equation Painleve-2 with zero 
in the right hand side of the equation (||) was investigated in works @, §|. In this case 
a solution in an interior layer near a bifurcation value of the parameter t* is determined 
by the equation Painleve-2, but already without a small parameter; and the problem, 
generally speaking, does not become simpler. 
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The asymptotics of the solutions for the Painleve equations with a leading term as an 
elliptic function with the modulated parameters were studied, for example, in works || 
pCfl . We must mention some of the works about the scaling limits or double asymptotics 
for the Painleve-2 equation. The general approach to the scaling limits of the Painleve 



equations based on the Backlund transformations was studied in [21]. The scaling limit 



passage from the equation Painleve-2 to Painleve- 1 was studied on the level of classical 



solutions in the work [22]. In work [£3| different approaches to the double asymptotics were 
developed. The qualitative analysis for the relation of the algebraic and fast oscillating 
asymptotic solutions of the equation (||) was done in the work [24] also. However, the 



asymptotic solutions constructed by this way are non-uniform with respect to two variable 
t and e. 

The major difference of the presented work in comparsion with the others cited is con- 
structing uniform asymptotic solution with respect to two parameters t and e as e — > 
into an interval of t, where the main term of asymptotics (elliptic function) is degenerated. 
This uniform asymptotic solution is valid on a segment containing the saddle-center bi- 
furcation point. In this work we specify a different types of asymptotic approximations of 
the being studied solution, their valid intervals and orders of neglected terms which result 
when these asymptotics are substituted into a being solved equation. 

The contents of the various sections are as follows. Section 2 states the main problem. 
The results are formulated in Section 3. An asymptotic expansion to be valid before 
the bifurcation point t* is written in Section 4. Section 5 is devoted to inner asymptotic 
expansions which are matched with each other and with asymptotics from Section 4. A fast 
oscillating asymptotic expansion of Kuzmak-type, a degeneration of the oscillations and 
the matching of the fast oscillations with inner expansions from Section 5 are described in 
Section 6. Open problems are discussed in Section 7. 

2 Naive statement of the problem 

Let's consider a cubic equation: 

2u 3 +tu = l, (2) 

which is obtained as a rejection of the term with the small parameter in the equation ([!]). 
There exist the point t* and the value u* such, that if t = t*, then u* is double root of the 
equation (^). The values u* and i* are easy to obtain by solving of the equations: 

2ul + Uu* = 1, Qui + U = 0. 

There are t, = -3 2^ 1 / 3 , u* = -4" 1 / 3 . 

The discriminant of the equation (|2|) has the form: 

t\ 3 fv> 2 



The discriminant D < when t < t* and hence the cubic equation (||) has three real 
roots tti(t) < U2(t) < Us(t). If t > then D > and the cubic equation (|2|) has one 
real root and two complex conjugate roots. At t = t* the roots u\{t) and 7/2 (t) coalesce 
Ui(i*) = u 2 (i*) = u*. 
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When t < it is possible to construct a real formal solution of the equation (|l]): 

oo 

u{t,e)=Y,e 2k u{t) (3) 

k=0 



, 



by taking any of roots Uj(t) as a leading term u(t). These three formal solutions are slowly 
varying equilibrium positions for the equation (|]). 
Consider an equation 

eV + (6u 2 j (t) + t)v = 0, j = 1,2,3, (4) 

which describe small perturbations of the leading term Uj(t) of the asymptotic expan- 
sions (||). If j = 2, then Qu^t) + t < 0, so that the equation (|j) has one exponen- 
tially growing solution and one exponentially decreasing solution. Hence corresponding 
asymptotic expansion is unstable with respect to small perturbation of the leading term. 
Otherwise if j = 1, 3 then the coefficients in the equation (Q) are positive and the formal 
solution (|3|) with Uj(t) as the leading term of the asymptotics is stable. 

When t = t* the roots u\(t) and U2(t) coalesce: = U2(t*) = and when 

t > i*, there exists only one equilibrium u-s(t). Therefore t* is the bifurcation point for 

the asymptotic solution (|3|) in case the leading term u(t) = ui(t). 

Our proposal is to construct a smooth asymptotic solution of the equation (|l|) with the 
leading term u\(t) when t < on a segment [t* — a, t* + a], a = const > 0. 



3 The main results 

Here we describe a smooth asymptotic solution constructed in this work. Following by 
V P Maslov |25j we will use the words "asymptotic solution with respect to mod (O (e a ))" , 
namely, a function is said to be an asymptotic solution of mod (O (e a )) of the equation ([!]) 
if after its substitution into this equation the latter is satisfied up to the terms of the 
order O (e a ). 

When t* — a < t < t*, (a = const > 0) and (i* — i)e -4 / 5 3> 1 the asymptotic solution 
with respect to mod (o (e 6 ) + O (e 6 (t - t*)" 13/2 ^ has the form 

u(t, e) = Ul (t) + b 2 "ffi® + e* 2 u(t). (5) 
The last term of the formal asymptotic solution as t — > t* — can be written as: 

lit) = o([t-uy 9 ' 2 ) . 

When \t — 1*| <C 1 the asymptotic solution is defined by two different types of the formal 
asymptotic expansions. The first one has the form 



u (t,e) =u, +£ 2 / 5 V)+£ 4 / 5 i(r). 



(6) 
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Here the variable r is denned by the formula r = (t — t*)e 4 / 5 , the function v(t) is defined 
as the solution of the equation Painleve-1: 

d 2 v(r) o 2 

— j-= 1- 6u*v +ti*r = 0, 

dr A 

with the pure algebraic asymptotic behavior as r — ► — oo: 



8w 



The formula @ is asymptotic solution with respect to mod (O (e 8 / 5 r 2 ) + 0(e 8 / 5 )) as 
K -r < e~ 4 / 5 . 

The function v(t) has poles of second order at some points r^, k = 1,2,... (see, 
e.g. @|): 

Near the poles the last term of the asymptotic solution can be written as 
v(t) = O ((r - T fc )" 4 ) , as r -> T fc . 

The expansion (||) is suitable at e _1//5 |r — 3> 1. The formula (JgJ) is asymptotic solution 
with respect to mod (o (e 8 / 5 ) + O (e 8 / 5 r 2 ) + O (e 8/5 r fc (r - r fe )~ 8 ) ) . 
As r — ► oo the main term of the asymptotics (y) is 

v(r,e) = Vrp{s, g 2 ,gs) +0 (t~ 7 ) , 7 = const > 0. (7) 

Here 

s = -r^ + Q a( X ), where x = e 2/5 ^ 7 / 4 , 
5 7 

the phase shift <r is defined in Section 5.4.1 by a formula (^). The parameter of the 
Weierstrass elliptic function g 2 = — 2u*. The second parameter #3 is defined by a solution 
of an equation (see, e.g. pT|): 

Re / dX uj = 0, 



where T is any circle on an algebraic curve uj 2 = A 3 + A/2 — 53/4. The last term of the 
asymptotics (||) has the form: 

v(t) = 0(t) + O [ -7 ) as r — ► 00 and t 7^ 17.. 

\{r-T k ) A J 



The formula @ is asymptotic solution with respect to mod^O (e 8 / 5 ) + O (e 8 / 5 r 3 / 2 ) + 
( e8/5 ' r fc /2 ( r -' 7 *) -8 ))- The expansion (|) is suitable as r « e" 4 / 5 and e" 1 / 5 ^ - 



rA^ !l » 1. 
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The second one, which is valid in the neighborhoods |r — (t^I 1 / 5 <C 1 of poles Tk of 
the function v(r), reads as 

u(t,e) = u* + w(9 k )+e 4 / 5 w(9 k ), (8) 

where 9 k = (r — T k )e~ 1 ^ 5 + e 1 I 9 k . The phase shift 9 k is defined by a formula (^)- The 
function w(9 k ) is defined by the formula: 

w{9 k ) 



4 + l6u 2 Jl 

The last term of the formal asymptotics (|8|) at \9\ — > oo can be written as 

i(0 fc ) = o(^N). 

The formula (||) is asymptotic solution with respect to mod (O (e 8 / 5 ) + O (s 8 ^ 5 ^^!)) • 

When (t — t*)e _2//3 3> 1 and i < t* + a the asymptotic solution with respect to 
mod (O (e 2 ) + O (e 2 (t — t*)~ 3 )) has the fast oscillating behavior: 

o 1 

«(t,e) = C/'(ii,t)+eC/'(ti,t). (9) 
Here the last term of the asymptotics (Mj at t — > i* + can be written as 

Z7(*l,f)=o((t-t*)- 3/2 ). 
The leading term of the asymptotic solution satisfies the Cauchy problem: 

(S'f [d tl U) =-U 4 -tU 2 + 2U + E(t), U\ tl=0 = u*. 



Here t\ = S(t)/e + 4>{t). The function E(t) is defined by the equation 

»a(t) 

'/3(t) 



/■<*(*) y 

I = 2 / ./-x 4 - ix 2 + 2x + dx = 2vr, 



where a(t) and (a(i) > (5{t)) are two real roots of the equation — x 4 — tx 2 + 2x+E{t) 
0, other roots of this equation are complex. 

The phase function S(t) is the solution of the Cauchy problem: 

_ /■"(*) A r 
T = S'V2 - S\ t =u =0. 

7/3(t) V-x 4 - te 2 + 2x + S(t) 

Where T is the constant defined by the formula 



2|uJV2 16-2P 
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where k ~ 0.463 is the unique solution of the equation 
-ky + k 2 + l 5/2 



o 



dy — xioU = o, 



and 

dy 



C,(fc) 



Vy[(y-fc) 2 + i] 

The phase shift 4>(t) is defined by an equation (see |2 

— — 4> = a = const. 
d E S' 

Remark 1. In this work the constant a remains out side of our analysis. Its value may 
be defined by using the monodromy-preserve method for the Painleve-2 equation |28[| . 

Remark 2. The domains of validity of the asymptotic solution (||) and the asymptotic 
solution (||) intersect, so that these expansions match. The solution of the Painleve-1 equa- 
tion which defines the asymptotics (g) has infinite sequence of the poles Tk, k = 1,2,.... 
Near all of these poles we match the asymptotic solutions (|6|) and (|J). As the number of 
the pole k — > oo the domain of validity of this complicated combine asymptotics (||) and (|8|) 
intersects with the domain of validity for the fast oscillating asymptotic solution (|9|). Its 
allows to match the sandwiched asymptotics with the fast oscillating asymptotic solution. 

4 The outer algebraic asymptotics 

The algebraic asymptotic solution (|5[) of the equation ([l]) is constructed here. This asymp- 
totic solution is suitable when t < t* and asymptotic behavior of this solution is investi- 
gated as t — > — 0. 

We construct the asymptotic solution of the equation (||) as: 

u(t, e) = u(t) + e 2 u(t) + e 4 u(t) H . (10) 

Let's formulate the result of this section. The asymptotic solution (|5|) with respect to 
mod (O (e 6 (t — i*) -13 / 2 )), where u(t) = ui(t) is least of the solutions of the equation @, 
is suitable when (t* — i)e -4 / 5 3> 1 and t > t* — a, where a = const > 0. 

4.1 Constructing the algebraic asymptotic solution 

Let's obtain the coefficients of the asymptotics (|To|). Substituting the ansatz (10) into the 



equation (|l|) and equating coefficients at identical powers of e we find the sequence of the 

k 

las for u(t), k = 0, 1, 2, . . . . 

2u 3 (t) + tu(t) = 1, (6n 2 (t) + tj u(t) = -u"(t), 
(6n 2 (t) + tj u(t) = -6u(t)i 2 (t) - u"(t). 
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The cubic equation for u(t) when t < t* has three real roots u\{t) < u 2 (t) < u^(t). As 



the leading term of asymptotic expansion (10) we choose u±(t). The second derivative of 
u(t) has the form: 

/ o \ ' o 
o„ _ ( u \ 2tu 2 (t) 



it 



3 ' 



1 



This allows to obtain the formula for u(t). 

It is easy to get the expressions for the following terms of the asymptotic solution (|l0"|). 
In an explicit form they are not adduced here, however, it is important to note, that the 

power of the denominator (^6u 2 (t) + tj in the coefficients of the asymptotics grows with 

each next step. The n-th term of the asymptotic expansion as (^6u 2 (t) + tj — > has the 
form 

u(t)=o((6u 2 (t)+t)~ 5n+1 ) . (11) 



Let's write the asymptotic behavior of the asymptotic expansion fllO| ) as t — > t*. For 
this purpose we shall calculate the asymptotics of the expression (^6u 2 (t) + tj : 

,0 2/ . 1 A n: r. r 2,. , N 5 



6u 2 (t) +t) = -2u*v / 6Vi7~ r I+ -(** -t) ^{t - U) 3/2 + O ((** - t) 2 ) . 



Using this formula and u(t), u(t) we obtain: 

u(t,e) =u„- —= Vt* - i + — ^— (t» -t)+e 2 
V6 18u* 



^(^-r 2 -o((**-r 3/2 ) 



+ o(e\u-t)- 9 ' 2 )+o({U-t) 



n3/2 



3 2 10 /3 



4.2 The domain of validity of the algebraic asymptotic solution 

The domain of validity for this expansion as t — > i* — is determined from the relation 
e 2 "?/ 1 (t) /u(t) <C 1. It follows from the formula (|ll]), that the expansion (|i~0| ) is suitable 



when (t* - t)e" 4 / 5 > 1. 

Evaluate the residual which is obtained when one substitutes the asymptotic solu- 
tion (||) into the equation ([[]) 

iT(t, e ) = _ e 6 C 2» + 2 1 2 + 6 1 3 \ _ £ 8 / g 2 2 + 1 2 2X _ ^£3 



fc 

Using the asymptotic behaviour of the u, k = 0, 1, 2 as t — ► — one can obtain 



F(t,e) = o(e 6 (t-t*)" 



13/2 
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5 The inner asymptotics 

In this section the asymptotic expansions of solution of (Q) which are suitable in the small 
neighborhood of a point i* are constructed. By following terminology of the matching 



method [29|, they are called "the inner asymptotic expansions". 



5.1 First inner expansion 

It follows from the consideration of the validity of the outer expansion, made in the 
previous section, that it is natural to make the following scaling of variables: 

(u - = e 2 / 5 v, (t - U) = e 4/5 r. 

As a result we write the equation ([!]) as 

^ + 6u*v 2 + u*t = -e 2/5 (to + 2v 3 ) . (12) 

In the limit as e — > we obtain the equation Painleve-1. This asymptotic reduction is 
known as one of the scaling limits for the Painleve-2 equation |3(| (see also [||, |l9|). 
A solution of this equation has the asymptotic expansion as r — > — oo: 

1 49 



v(r,e) = ( -^r/6 + ^ + 768 ^ ( _ r)9/2 + 

+ e 2/5 ( _ J_ + 1 +...Uof e^r^ 2 



18u* 144 v / 6(-r)3/ 2 
The asymptotic solution of the equation (1121) we build as: 



v(t,e) = "v(t) + V^ 5 S(r), (13) 



n=l 



where the function v(t) is the solution of the Painleve-1 equation. 



Here it is shown, that the asymptotic solution ( |I3| ) is suitable in the neighborhood of 



infinity (when — r <C £ 4//5 ) and in the neighborhood of the poles for the function v(t): 
(r-r fc )e- 1 /5 >>1 . 

5.1.1 Asymptotic behaviour as r — > — oo 

The coefficients of the asymptotics are calculated from the matching condition for the 



asymptotic expansion (10) as t — > t* and the expansion (13) as r — > — oo. In particular, 
v(t) has the algebraic asymptotics: 

g (T )| T _ = -^ +i J-, + ^ + 0(r-'). (14) 



In the book [26| it is shown, that there exists the solution of the Painleve-1 equation 



with the asymptotics (uA)- The data of a monodromy for the solution of the Painleve-1 



equation with the asymptotics (14) are calculated in the work 



74 



O M Kiselev 



The first correction in the asymptotics (|13| ) satisfies the equation 

»i 

'v 

dr* 

The asymptotics of the solution for this equation as r — ► — oo has the form 



d 2 v 01 o „0q 

+ 12u*vv = -tv - 2v 6 . (15) 



v(t) = — + ^ + O (r" 4 ) . 

V ; 18u* 144^6(-r)3/ 2 V ; 

Asymptotics of the higher corrections is constructed by ordinary way. The n-th cor- 
rection as t — ► — oo has an order: 



S(r) = O ((-r)(" +1 )/ 2 ) . 



5.1.2 Validity of the asymptotic solution as r — ► — oo 

The requirement of validity for the asymptotics is e 2 ^ 5 v/v <C 1. It reduces to the condition 
(-r) « e~ 4 / 5 . 

The residual of the asymptotic solution has the form: 

F ( T) e ) = _ e 8/5 + T l + gl ^ _ £ 10/ 5 1 2 _ £ 12/5 1 3 _ 

Using the asymptotic behaviour of v, k = 0, 1 as (— r) <C e _4//5 one can obtain 
F(r,e) = 0( £ 8 / 5 (r 2 ^ 



5.1.3 Asymptotic behaviour near the poles 

The function v(t) has the poles when r G (—00,00). Let's denote these poles by r^. In 

the neighborhood of the pole r — ► ± the function S(r) is defined by the converging 
power series (see e.g. 

(16) 

The constants and are the parameters of this solution. In the review [11] it is noted, 
that the problem on the connection between the asymptotics of this solution at infinity and 
the constants Tk and is not investigated yet. The points of the poles and appropriate 
constants can be obtained with the help of the numerical calculation using the given 
asymptotics at infinity (|l4|). 

The asymptotics of v as r — > =t may be written as a sum of a certain solution of a 

nonhomogeneous linearized Painleve-1 equation v c (t) 

*M - -J^ ^~<T- r k) + - r k? + 0((r- r k f) , 
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and two solutions of a homogeneous linearized equation ui(r), V2(t): 

Mr) = 7 ^ + - r k ) + ^(r - T k f + 2c k u*{r - r fc ) 3 + O ((r - r fc ) 5 ) , 

(t - T k y 10 5 



^(r) = (r-T fc ) 4 + 0((r-T fc ) 8 ). 

Thus: 

i = i c (r) + a^i(r)+6^ 2 (r). (17) 

i, i. 
Here a fc and 6 fc are constants. 

Higher corrections have the same form: 



u(t) = u c (t) + a fc «i(r) + b k v 2 (r) 
where 

^ c (r) = 0((r-r fc )- 2(n+1) ), 



as r — > Tfc. 



5.1.4 Validity of the asymptotic solution as r — ► Tfc 



By using the asymptotics fll6| ) and fll7| ) it is easy to see, that the asymptotic expansion (13) 
is suitable at 

e- 1/5 \r-T k \ » 1. 

The residual of the asymptotic solution as r — > when e~ 1 ^\r — T k \ S> 1 is 



^(r, e) = e 8/5 ( (^4) + ((r - r,)" 8 ) . 



5.2 Second inner expansion 



For the construction of the uniform asymptotics in the neighborhood of the pole of the 

function v it is ne 
function (see ||): 



function v it is necessary to make one more scaling of the independent variable and the 



-2/5 



e 1 v = w. 



For function w we obtain the equation: 

^ + 6u*w 2 + 2w 3 = -e 4/5 T fc (u* +w)- e9(u* + w). (18) 
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The solution of this equation has following asymptotic expansion as 9 — > — oo: 

/i- \ 



1 1 

w = — 



„^ + 4# + °( e - 8 ) +£l/5 (w + ( 9 " 4 ) 

1_\ 2 1 



U, 



:9 4 



+ £4/5 /_ 12O^ + ^ 2 + o(0 -i^ 



(19) 



1. 



2 



+ e 6 / 5 ^0 2 + c fc fl 4 + ^0 2 + O (e 1 ) ) + O (e 7 N* 



+ e 



8/6 6 



1 T 2 7 ,3 
4 T fc "* />6 



300 



6 + o(0 2 )) + o(>/ 5 ). 



l 

This long asymptotic formula shows, that the constant appears only in the correction 
of an order e 8 / 5 . If we want to construct the first correction of the asymptotics for the 
first inner expansion after the pole r k , we must construct the correction in order e 8//5 for 
the second inner expansion. 

It is convenient to include a time shift depended on e into the main term, and construct 
the asymptotic expansion depended on a new time variable: 

n 

where 9 k = const. 

We search the asymptotic expansion for the solution of this equation as a segment of 
an asymptotic series 

w(9 k , e) = w(9 k ) + e^w^k) + ew(9 k ) + e^w(9 k ) + e^w{9 k ). (20) 
In this case the equation for the w(9 k ,e) looks like: 

6u*w 2 + 2w 3 = -e 4 ^T k {u* + w) — e9 k (u* + w) 



d 2 w 



del 



+ e 6/5 9 k (u* + w) + e 8/5 9 k (u, + w) + 



It is shown here, that the asymptotic solution (0) is the formal asymptotic solution 
of the equation |[|) with respect to mod (O (e 8/5 r|6» 4 ) + O {e 9 / 5 T k 9 5 ) + O {e 2 9 6 )) when 
\9rl /5 \ «e 1 /5. 
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The solution of the equation for the leading term of the asymptotics ( j20[) is defined by 



the asymptotics as r — > of the asymptotic expansion (13), which is outer with respect 
to (pO). This solution has the form 



16u* 



4 + lbu 2 8 2 ' 



(21) 



The constants 9 k are defined by asymptotics of the function w(9,e). Using the for- 
mula (Ii~9|) we obtain: 



n= 1,2, 



The corrections in the expansion (120) satisfy the linearized equations 



(22) 



d 2 w 



Jo\ 1 



+ \2u*w + bw ) w = —r k lu* + w), 



d w 



+ (l2u*w + 6u; 2 ^ w = 9 k (u* + wj , 



d w 

~dW 



^0o\3 \ ( 0\ 
\2u*w + bw w = 9 k lu*+w , 



d 2 w 
Iff 2 



^0 2 \ 3 
+ \2u*w + bw w 



bw 2 [w + tt* ) + 9 k [u* + w 



The expression for w can be used to obtain two linearly independent solutions of the 
homogeneous equation for the corrections: 



»-'i 



2 ' 



W 2 



+ 2ulel-uji + -ei + 



2ui 
7 



(1 + 4^) 



2 ' 



By using these solutions of the homogeneous equation it is easy to get the solutions of 
the nonhomogeneous equations for the corrections. The asymptotics of the corrections as 
9 — > oo has the form: 



i = + _Jk_ + n_ e - 2 



1 



10 



120u* 120 



rf 4 + ote«). 
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3 

It is important to note that the leading term of the asymptotics of w as 9k — > ±oo 
is the same. This term defines the constants Ck and hence the solution of the Painleve-1 
equation before and after the pole. We have the same value of the constant c& as 9k — > — oo 
and 9k — > oo and then we have the same asymptotic solution of Painleve-1 equation before 
and after the pole 

4 

An opposite result takes place for the coefficient w in order 9\ as 9k — ► oo. This 

i_ i_ 
coefficient is changed. It is equal b k / (56it*) as 6>fc — ► — oo and 6 fe / (56ii*) — lltt*rf/2100 

as 0fc — > oo. 



5.2.1 Validity of second internal asymptotic solution 

Using the asymptotics for the corrections and the leading term, we obtain, that the ex- 
pansion ( |2C)| ) is suitable when \9t^ 5 \ <C e -1 / 5 . On the other hand, the expansion (13) 
is suitable at \9\ 3> 1. Hence, the domains of the applicability for the expansions (13) 
and pi) are intersected at realization of the condition | | e . If we take into account 



also the requirement of fitness of the asymptotic expansion (20), then get the restriction 
\T~k\ "C e _4//5 . From this inequality it follows, that in this section the formal asymptotic 
expansions to be suitable when \t — 1*\ -C 1 are constructed. 

We calculate the residual of second internal asymptotic expansion using the asymptotic 

k 

behaviour of w, k = 0, 1, 2, 3, 4: 

F(9,T k ,e)=e 9 / 5 0(T k 9 5 ), (23) 
when I^Hrfcl 1 / 5 < e" 1 / 5 . 



5.3 Dynamics in the internal layer 

Using the asymptotic expansion of the second inner expansion as 9 — > oo and the first 
inner expansion as r — ► r k + we find that the first inner expansion after the pole has the 
form ( |l~3| ) where 



v(t) = v c (t) + a£vi(r) + b^v 2 (T). 

Here 

n , n _ n n_ n 

< = a fc< K = b k+ A k, n=l,2. 



The shift may be calculated from the asymptotics of the second inner expansion as 
9 — > oo, for n = 1 we have obtained: 



75 

Thus a behaviour of the asymptotic solution in the internal layer is combined by the 
first and the second inner asymptotic expansions. 
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5.4 The asymptotics of the inner expansions as r — > oo 

In the above sections we demonstrate the asymptotic behaviour of the asymptotic solution 
at t — ► — co and near the poles of the solution for the Painleve-1 equation. Below we study 
the asymptotic behaviour of the solution as r — > co. 

The regular asymptotic expansion on e is constructed in the previous section concerning 
the first inner asymptotic expansion. This asymptotics is not valid as large r. To use the 
first inner asymptotic expansion as r — ► oo we must modulate the parameters of the 
solution of the Painleve-1 equation and cancel secular terms in the first and the second 
corrections in the asymptotic expansion (|j~3|) using singular perturbation theory. Instead of 
modulating the parameters of the solution of the Painleve-1 equation it is more convenient 
to study the modulation equation for parameters of the main term of the asymptotics for 
the solution of the Painleve-1 equation as r — ► oo. This study will be developed in this 
subsection. 



5.4.1 Asymptotic behaviour of first inner expansion 

The elliptic asymptotics of the solution for the Painleve-1 equation as r — ► oo was obtained 



by P Boutroux [12|. Here we are interested in a connection formula for the solution of the 
Painleve-1 equation. Namely we have the asymptotic behaviour of the solution as r — > — oo 
and we need the asymptotic behaviour of the same solution as r — > oo. The Painleve- 
1 equation is integrable by the monodromy-preserving method |28|j . If the monodromy 
data are known, then the solution of the Painleve-1 equation is uniquely defined. In the 
correspondence with ||], in our case the monodromy data are constants S2 and S3 and 

these constants are equal to zen 
poles has the form (see e.g. [|l]| ) 



these constants are equal to zero. The asymptotics of the function v(t) outside of the 



8 = ^ 3( ff ) + O (r-T) , (24) 

where 7 > is some constant, the function p(cr) is determined by the Weierstrass elliptic 
function 

P( a ) = ~P( a i 92,93 )/u*- (25) 

The phase function a = |t" 5//4 - It is important to note, that in the formula (|24| ) the shift 
of the phase function a is equal to zero. A parameter is 52 = — 2n^ and a parameter 53 is 
defined as a solution of the equation: 

Re / u dX = 0, 

where 7 is any circle on an algebraic curve: to 2 = A 3 + A/2 — 53/4. 

We are interested in the asymptotic solution of the perturbed Painleve-1 equation (|i"2|). 
Therefore the asymptotics of the works ]i~9| , 11 1 for unperturbed Painleve-1 equation are 
the main term of our asymptotics on e. Substitute into (|l~2^): 



v(t,s) = y/rp(a,e), 



(26) 
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where a = 4r 5 / 4 /5. As a result we obtain an equation: 



p" + 6u*p 2 + u* 



5a 

T 



5a Y 2 



2/5 



5(7 



2/5 



(p + 2p J 



(27) 



Let us to construct the asymptotics of the perturbed equation as a segment of the asymp- 
totic series: 



P (a,e)=p(s)+e^ (|) 



2/5 



l p{s)+e^{°- 



K \ 4 /5 
00" \ ' 2 



(28) 



where 



s = a + a(x), 



2/5 5 /5a \ 7/5 



7 v 4 y 



o . 



a is modulated phase shift and x is once more slow variable. 

Substitute the formula (|28| ) into the equation (p7|). Then let's equate the coefficients 
with identical powers of e. As a result we obtain the sequence of the equations: 



U„ Uo 

p + ou*p + u* 



X o, 1/5 V'o 



(29) 



01 



p + Ylu^pp 



^o s \ 0,0„ 
p + 2p 3 -a'p" 



+ 



1 /5 



2 V 4 



0/0, 

o-p 



i, 7 /5 
P + < 



(30) 



p + Xlu^pp = — 6u*p — p I 1 + bp 



o„ 0„o' 0,1,, 
p — a p — a p 



I £l (p,P,P 



4cr 1 L 2 (p,p,p 



(31) 



Here Li 2 are linear operators. 

Solutions of these equations may be represented as asymptotic expansions as a — > 00. 

We will assume that the corrections of these expansions are small as a — > 00. Therefore 

we will neglect the correction terms of this asymptotics and consider only the main terms 

012 

of these asymptotics with respect to a. The equations for the main terms p , p , Pq have 
the form: 



p + 6u*p + = 0, 



(32) 



Po + 12u *PoPo 



(V ^0q 

(Po + 2po 



o,o„ 
° Po, 



(33) 
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2// , 10 2 1 2 1 / 2 \ f0,\ 2 0„ 0„0, 0,h, 

p + 12u*p Po = -6-u*p - p (1 + 6p J - [a J p - a p - a p (34) 

The requirement of validity for the first inner asymptotics at r — ► oo is e 2n / 5 o" 2n / 5 p/p <C 

1, where n = 1,2. We find the modulated equation for the a such that the first inner 
expansion is valid as large r. 

The solution of the equation for the main term is the function p(s) which is defined 
by (|25|). The equations for the corrections is the Lame equations with external force. To 
write the solutions of these equations we will use two linear independent solutions of the 
Lame equation. Denote one of these solutions by 

pi(s) = d s p . 

We denote the second solution as P2(s)- The solutions pi(s) and P2(s) are such that a 
wronskian: 

W(p 1 ,p 2 ) = 1. 

The second solution is aperiodic: 

p 2 (s + Q) = Cpi(s) +p2(s), where C = const ^ 0. 

The solution of the equation ( |33[ ) can be written as: 

Po( s ) = A k p!(s) + B k p 2 (s) +pi(s) j dz (-p (z) - 2p%(zf) p 2 (z) 

(35) 



p 2 (s) / dz(-p (z)-2pl(z)jpi(z)-a'sp 1 (s). 

JSQ V 7 



Here sq = Sk + 0/2, where s k is pole of the function Pq(s) and $7 is a real period of the p. 

l l 
The A k and B k are constants. 

The first correction p is bounded as s £ R if 

k = ^' + ^ R-P- [ dz (J> (z) + 2°pl(z)) Pl (z). 
In this formula the integral must be regularized. Namely: 

i?.P. y dz (p (z) + 2/>g(2)) pi(z) = res r=0 - J dz (p (z) + 2?§(z)) p x (« 



Using a perturbation theory for the second order equation developed in works [32| — [27] , 

2 

we can show, that the solution for the p has an order 0(a), a ^ <j k , k € Z, if the function a 



o, 



l 

is a solution of the differential equation: a" = 0. This equation and the equation for the B k 
allow to write the Cauchy problem for the a in the form: 



i 

a' = CB k + R.P. dz (p (z) + 2p 3 (z)) Pa(4 



(36) 



82 



O M Kiselev 



where 

Xk ~ £ — ) ■ 

1 1 

The constants and are defined by matching conditions of the asymptotics (|26|), 



H) and dH) as r fe -> r fe + 0: 

A fc = ^ - A.p. r /2 ^ (? (z) + 2?§(z)) P2 (z), 

n Jo v y 



5* = ^ - dz (? (z) + 2?g(z)) P i(z). 

As a result we obtain the formula (0) for the main term of the asymptotics 

5.4.2 Validity of the first inner expansion as r — ► oo 

Let us denote by 

« (r) = T (n + 1)/2" ((T); n = ,l,2. 

Then the condition of the validity e 2 ^v/v <C 1, n = 1, 2 fulfills as e 2//5 r/\/r C 1 or as the 
same: 

T« £ - 4 / 5 . 

Near the pole r n we obtain the asymptotics: 
g = O ( , , and h = 



T ~ T k) 2 J ' \( T - T fc) 4 

Hence the first inner asymptotic expansion is suitable near the poles as 

e - 1/5 (T-T fc K: 1/4 »i. 

The residual of the first inner expansion is 

F(t, e) = —e 2 (l2u*vv — 6vv 2 + 6v 2 v + ri?j — e 12 / 5 (^2v 3 + 12vvv + 6u*v 2 
+ e 14 / 5 (d 2 l + 6°' 2 ) + e 16 / 5 6^ 2 + e^2p. 

Using the results of the above section outside of the poles of v one can obtain: 



F = e 2 r^ 2 )+0 



2,5/2 



£ T 



10 
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5.4.3 Validity of the second inner expansion as — > oo 

In the second inner expansion the asymptotics as r — > oo corresponds to the asymptotics 
at Tfc — > oo. It is easy to see, that in this case the first correction grows. This growth 
limits the value Tfe, at which the asymptotics (|20|) is correct: |t&| <C e -4 ^ 5 . One can use 



the formula (23) to obtain the residual of the second inner asymptotics as — > oo. 



6 Fast oscillating asymptotics 
6.1 The Kuzmak's approximation 

In this section we apply formulas obtained in |H]]-[j34|], |^7j to the fast oscillating formal 
asymptotic solution of the Painleve-2 equation. These formulas are usable when t > t*. 
The fast oscillating asymptotics is constructed as 

o 1 2 

u(t, e) = U(ti,t) + eU(t u t) + e 2 U(h,t) + ■■■ . (37) 

As the argument t\ we use expression 

h = S{t)/e + (j){t), 

where S(t) and 4>(t) are unknown functions. 

The equations for the leading term and the corrections of the asymptotics ([}?]) look 
like: 

(S'fdlU + 2U 3 + Ut = 1, (38) 



(S'fdl U + ( QU 2 + t)ll = -2S'd 2 ttl U - S"d tl U - 2S'(j)'dl U, (39) 



{S') 2 d 2 Jj+ (qU 2 + t^jU = -6UU 2 - 2S'd 2 ttl U- S"d h U - 2S'(j) , d 2 1 U 

- d^U - ((p') 2 d tl U- <j)"d tl U - 2<f>'d t d tl U. 
o 

Integrating once with respect to t\ the equation for U we obtain: 

„ / o\ 2 , „ 
{S') 2 (d tl U) =-U i -tU 2 + 2U + E{t), (41) 



(40) 



where E(t) is the "constant of integration". 

We study the equation when the right-hand side has two real roots 0(t) < a(t). We 
define the initial data as 

U\ti=o = Pit)- 
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Solution of the equation (|4l| ) is an elliptic function. The left hand side of the formula (41) 

o 

is positive for real functions S' and U. The term of highest order in the right hand side 
o 

is — U 4 . Therefore the real solution of the equation (|4l| ) has no poles when t\ £ M and 
\E(t)\ < oo. 

To construct an uniform asymptotics (37) we choose the unknown functions S(t), <f>(t) 

and the "constant of integration" E(t) by a special way. They must satisfy to anti-resonant 

conditions for equations (|38|)-(^0|). The first condition is the boundedness of the right hand 

o 

side of ( p9| ) as t\ G E. It is satisfied if a period of the oscillates of the function U{t\,t) 
on t\ is a constant (see for example ISjJ ): 



T = v 7 ^' / ; (42) 

J/3(t) \/-x 4 - tx 2 + 2x + E(t) 

l 

The next condition is a boundedness of the first correction U(t±,t) when t\ G R. It 
gives equation defining a main term of an action Iq (see, |3jJ ): 



h = S' 



T 



dtMti,t) 



dt\ = const. 



Using the explicit expression for the derivative with respect to t\ we present this formula 
in a some other form: 

a(t) 



Jo = 2/ V-x 4 - tx 2 + 2x + E(t) dx = const, (43) 

J/3(t) 

where a(t) and (3(t) are the solutions of the equation — x 4 — tx 2 + 2x + E(t) = 0. The 
equation ( f43| ) means that the action is a constant. At the other hand-side the equation (43) 
defines the energy E(t) as a function with respect to slow time t. 

2 

At last the necessary condition of boundedness of the second correction U with respect 
to t\ is the equation (see, 



M^' = a = const. (44) 

The equations (f42^)-(|44|) define the parameters of the main term of the asymptotics 
as functions with respect to t. To find these functions one should define corresponding 
constants. 

Notice that to define the value of the action Jo we construct the asymptotic solution of 

o 

the equation (||) when t > t*. Thus the polynomial of the fourth power on U in the right 
hand side of the equation ( |4"T| ) can have no more than two various real roots a(t) and /3(t). 
Hence this polynomial can be submitted as: 

F{x, t) = (a(t) -x)(x- /3(t)) ({x - m{t)f + n 2 (t)) . 

The degeneration of the elliptic integral at t = i* corresponds to the case m(t*) = 
P(t*) = u* and n(t*) = 0. For this case it is easy to calculate the constant in the right 
hand side of the equation ([43|), which is equal to 2ir (i.e. Jo = 2ir) and the value of the 
parameter E(t*) = E* = | (\) 2 ^ ■ 
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6.2 Degeneration of the fast oscillating asymptotics 

In this subsections we calculate the asymptotic behaviour of the phase functions S(t) as 
t -> U + o. 

The oscillating solution is degenerated as t — > t* + 0. Let's construct the asymptotics 
of this solution in the neighborhood of the degeneration point. For this purpose we cal- 
culate the asymptotics of the phase function S(t) and the function E(t). Let's write the 
equation (E3|) as: 



P 



y/(a-x)(x- 0) [(x - m) 2 + n 2 ] dx = ir, (45) 



where a, (3, m, n are real functions when t > i*. These functions satisfy the Vieta 
equations: 



a + P + 2m = 0, 
m 2 + n 2 + a/3 + 2m(a + 0) =t, 
(a + 0) (m 2 + n 2 ) + 2ma(3 = 2, 
af3 (m 2 + n 2 ) = -E. 



(46) 



The equation ((42|) and three equation from (|4q) define the dependency a, (3, m, n on 
the parameter t. The last equation in (^) defines the function E(t). Let's make changes 
of variables: E = E* + gi , t = t* + r?, 4m = + mi. After simple transformations of the 
equations ([46]) we obtain: 



2m* [6m? - 2n 2 + ?]] + [2m? - 2n 2 + 77] 2mi = 0, 

m 2 (12m? - 4n 2 + r/) + 2m^m 1 (6m? - 2n 2 + 77) (47) 
+ (3m? — n 2 + 77) (m? + n 2 ) = —51. 

Construct the solution of this system as t — ► + as: 

mi = ny/rj + 0(rj), n = i>iy/rj + 0(r/), 5l = 7l? 7 + O ( ?? 3/2 



Let's substitute these expressions in (|47|), equate the coefficients at the identical powers 
of rj. As a result we obtain: 

6/i?-2z^ 2 = -l, 7l =m 2 . 

To define the constants fii and ^1 it is necessary to construct the asymptotics as rj — > +0 
of the left hand side of the equation (p5|). The asymptotics of the outside the integral 
coefficient in the equation (45) has the form 



(a- (3y = 64|m 
The integral in the equation ([45|) is presented as 



(48) 



I{k,S)= / dz - z)zy/(z - k6) 2 + S 2 , 
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where 

x-P m-P 2 n 2 

a — p a — p (a — p)^ 

The value of the constant will be defined from an asymptotics below. 
The asymptotics of an integral I{k, 5) as 5 — > has the form 

J(*. 5) = ^ - M| + + C (fc)<* 5/2 + O (5 3 ) , (50) 

where 

' 5j "[fe-^ + l] 5 / 2 " 

First three terms in this formula are calculated by standard way. Let's show as we can 
obtain the function c(k). For this purpose the following trick (fl35[) is applicable. Let's 
calculate third derivative with respect to 5 of the function I(k,6): 

d 3 I „ f 1 , r- r kz + k 2 S + 5 



dP = ~ 3 Jo dZ ^^ X ~ z)z 



[(z-k6) 2 + 6] 5/2 ' 

On the right hand side we replace z by 5y and we present the integral as 



Solving the ordinary differential equation ( pl|) in the neighborhood of 5 = 0, we get: 

I(k, 5)=c + 6a + 6 2 c 2 + 6 5 / 2 ^-c 3 {k) + O {5 3 ) , 

15 



where 



oo 



Jo [{y-k) 2 + lf 2U 

After that it is easy to obtain the asymptotics (|50|). 

To define the value of k we substitute the asymptotics ( 48 ) and (|50D in (f45|) and equate 



to zero the coefficients at identical powers of rj. In the result we get at r/ 5 /^he equation 
c(fc) = 0. 

This is the transcendental equation for the definition of the parameter k. The numerical 



solution gives k ~ 0.463. Using the formula (49), we get: 



k 



Ml = 3 |iy| ' Ul = V W^~Y 
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To construct the asymptotics of S(t) as t — > t* + we use the equation connecting the 
period of fast oscillations with its phase (plf): 



T = V2S' dx (52) 

J 3 \ (a — x)(x — 3) \(x — m) 2 + n 2 ] 



1/3 \/ {a — x)(x — 3) [(x — m) 2 + n 2 ] 
Present the integral in the right hand side as 
If 1 dz 



a -(3 J — z)z \(z - k5) 2 + S 2 ] 

After the same replacements, as at the construction of the asymptotics as 5 — > we 
get: 

J= 6 ——\ dy = + Q(l). 

a-BJo ^fy [(y + fc)2 + 1] 

We substitute this expression into the equation ([52]) , use the asymptotics 8 and (a — (3) 
as r] — > +0 and in the result we get: 

S' = (t- U)V*S*(k) +of(t- U) 1 ^ , 

where 

* w = As* V /4 , aw ^ * 



V^2 V6-2A;V ' W 7o y/y [(y - k) 2 + 1] ' 

o 

The period of the oscillations for the function U(ti,t) with respect to the variable ii 
in the Krylov-Bogolubov's method is an arbitrary constant. Let's choose it such, that 
S*{k) = 1: 

g^)vga(fc) / 3 y /4 

r= 2Kp \j^2k 2 ) ■ (53) 
In the result the phase of the oscillations as t — > t* has a form 

S(t) = |(t - t,) 5 / 4 + O ((< - t*) 3/2 ) + So, (54) 

where 5o is some constant. Its value will be defined below at the matching of the asymp- 
totics (|37|) and inner asymptotics (fH|), (p0|) as i — > t* + 0. 

6.3 The domain of validity of the fast oscillating asymptotics 

In this subsection we establish the domain of validity of the fast oscillating asymptotics 
and compute the residual of this asymptotic solution. 

l o 

The validity of the asymptotics is defined by the formula eU U. Let us check this 
requirement. For this we must obtain the order of the first correction as t — > t* + 0. 
Evaluate the order of the right hand side of the equation for the first correction: 

F 1 (t 1 ,t,e) = -2S'dl 1 U-S"d tl U. 
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From the equation for U one can evaluate second term in Fi(t\,t, e) as t — > t* + 0: 

s"d tl u = o((t-u)- 1 ) . 

o 

One must reduce formula for the derivative of U with respect to t to evaluate of the 
first term in the formula for Fi(ti,t,e). 



The function U is the inverse function with respect to the elliptic integral 

dy 



h + 1 = S' I 
J a 



hit) V-y 4 - ty 2 + 2 ?/ + W) 

Both limits of the integration are functions with respect to t, it is not convenient for us. 
Make the substitution y = (a — (3)z + (3. Then we obtain 



h+to 



S' 





a-/3 



dz 



(a -13) J ^(1-^)^(^-7)2 + ^' 



Now we differentiate this formula with respect to t, as a result we obtain the formula for 

■ o 



the dt 



u-P 



dt 





u-p 

a — f3 



(a-(3)S'\ 



ON / 



a-U )[U-0 ) \ [U-m] +n 



,5" 



S' 





a — 13 



dz 



(a -[3) J Jo ^z{l-z){{z- 1 f + 5 2 ) 
dz 55' — (z — 7)7' 





u-P 

a-0 



(a-/3)7 y/z(l- z)({z-j) 2 + 5 2 ) 



3/2 



where 7 = 

' a— b 

In the same way we can obtain the formula: 



dt 





U-a 
a — (3 



1 



(a-/3)S'\ 



\ / 



a-U) [U-P] [ [U-m] +n 



-dt 



S' 



s 1 





U-a 
a-0 



dz 



(a-P)JJ_ 1 v / z (i_ z )(( z _r)2 + 52) 
dz 55' -iz- T)V 



(1 

U-a 



{a~P)Jo y/z(l- z)((z-T) 2 + 5 2 ) 



3/2 



where T 



a- /3 • 
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These formulas will be useful when we will reduce the formula for the second derivative 
o 

of U with respect to t. 

o 

The first derivative of U with respect to t has the form: 



o 



d t U =(a-P) 



a' - [3' 
(«-/3) 2 



+ d t 



-P 



a — (3 



+ 



(<x-0)S'\ 



\ / 



-dt 



S' 



S' 



a-U )[U-P ) \ (U-m) +n 



dz 



o 

u-P 

a-8 



(U-P)JJ ^(1 _ s) (( z _ 7 )2 + £2) 

dz (z - 7)7' + 55' 



o 

a — [3 



{a-P)Jo v / z(l-z)((z- 7) 2 + ^) 



Now we can evaluate the second derivative d^ ti U. 



3/2 



d 2 ttl u = d t 



i 

s 7 ^ 



o\ /o 



a-U] [U-P] [ (U-m) +n 



o 



Using the formula for dtU one can obtain as t — > t* + 0: 

dl 1 u = o((t-u)-^). 

This formula allows to evaluate the right hand side in the equation (^9|) as t — > + 0: 

F 1 (t 1 ,M) = o((t-g- 5 ^ 



The first correction is periodical function with respect to t%. One can derive the solution 
of the equation for the first correction using two linear independent solution of the equation 

{S'fd 2 tl V + [QU 2 + t^j V = 0. 

o 

Here our goal is to write these solutions in the terms of U because then we evaluate the 

order of derivatives of the first correction of asymptotic solution ([}?]) using the formula 
o 

for dtU. 

The first one is 



*7i(ii, t, e) = d tl JJ = ±lf y -m - til 2 + 2U + E(t). 



Here the sign before the root is +, when dt-JJ > and vice versa. 
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The second solution of the homogeneous linearized equation for the first correction is 

da 



1 



U 2 (h, t, e) = ±— Y ~U A - tU 2 + 2U + E(t) 



to 



-m - tU 2 + 2U + E(t) 



Integral in this formula must be regularized because integrand has second order poles 
oo 

at points U = a and U = (3. One of the possibile way of the regularization is done in |3q| . 
Here we will follow |36| . 

Near the singular points a = t\ j G Z one must represent the integral as 



S' 



da 



t j_ A Uf(a,t,e) 



-S' 
S' 



t{+\ 



1 



da 



t |_ A d a Ui(a,t,e) \Ut(a,t,£) 



1 



Ui(a,t,e)d a Ui(a,t,e) 



t{+\ 



S' 



da d 2 Ui(a,t,e) 



t{-\ U 1 (a,t,e)(d a U 1 (a,t,e)) 



2 ■ 



The parameter A may be for example A = T/4, where T is the period of oscillations of the 
o 

function U with respect to t±. 

Change the second derivative of U\ in the last integrand as 



i 



and the first derivative of U\ as 



1 



dtJJx = S'dtJJ = ^ I 1 - tU - 2U- 
As a result we obtain formula for regularization of the integral: 



S' 



da 



(S 



'\3 



t{-\ 



U 2 (a,t,e) 



o 



o 



Ut(a,t,e) ( l-tU(a,t,e) - 2U 3 (a,T,e 
o 



a=t{ +A 
cr=t[-\ 



+ s' 



t\+\ 



da 6U 2 (a,t,e) +t 



I l-tU(a,t,e) -2U 3 (a,t,e)\ 



2 ' 



Using the functions U\ and U2 one can solve the equation (39) for the first correction 



of the asymptotic solution ( |37| ) and obtain the solution the terms of U . 

1 

One can see the first correction U has the order of the right hand side of the equation 
for the first correction multiplyed on S'. It means that 



u = o ((t - ^r 3/2 

as t — > t* + 0. 
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This formula allows to obtain the restriction for the validity of the formal asymptotic 
solution @: 

{t-U)e~ 2/3 > 1. 

Evaluate the residual of the asymptotic solution (^). For this we must evaluate the 
function 

F(t lt t,e) = -e 2 d 2 U - e 2 d t ~ ^^t/ - e 3 d 2 U. 

For all t G + a] the order of F(ti,t,e) = 0(e 2 ). But the order of F grows as 

t — > i* + 0, because the derivatives with respect to t have singularity at t = t*. In the right 

o 

hand side of this formula there is only functions on the U, then we can differentiate the 

o 

right hand side for deriving of the second derivation of U with respect to t. The formula 

o 

for d 2 U will be very large if we will write it in these definitions but now one can evaluate 
o 

the order of d 2 U as t — ► £* + using the formulas 

S' = O ((t-t*) 1/4 ) , 0{ 1 ) = 0{5) = o{{t-U) 1/2 Y as t -> U + 0. 



As a result one obtain: 
o 



dfU(h,t,e) = 0({t-t 



-10/4 



Using the same formulas one can evaluate the order of the F(ti,t, e) as t — > t* + 0: 
F(t u t, e) = 0(e 2 (t- U)- 11 ' 4 ) +o(e\t- Q~ 17 ^ 



6.4 The matching of the fast oscillating asymptotic solution 
and the inner asymptotics 

The matching of this asymptotics with the inner asymptotics ( |l3| ) and (20) is carried 
out. From the matching condition for the phase function we obtain the initial condition 
S(t)\ t=tt =0. 

o 

Now we turn to the evaluation of the asymptotics for the function U as t — > i* + 0. The 
o 

function U may be written in the implicit form: 

h = -S' L - (55) 
Ju V-x 4 - tx 2 + 2x + E(t) 

We remind the parameter t\ is equal to e~ l S{t) + 4>(t) and the additional term So is 
undefined in the function S(t). This term we define in this subsection. 



The formula (55) allows us to obtain the main term of asymptotics for the U at t 
i* + 0. Denote 

o 

u = u* + w(h,t) 
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and 

x = -u* + y. 

Using asymptotics of E(t) and a(t) as t — > t* + one can get 

r«*+a(t) dy 



h = -s' / 

Jw 



\Jw V-y*-^*y 3 v ; v 



This formula allows to write the asymptotic expansion of U as t — > t* + in the form: 
t) = Wo(ti/5') + O ((t - t,) 1 / 2 ) + O (^^) ■ (56) 

The main term of the asymptotics is defined by formula: 

4u* 



l + 4n2(t 1 /5') 2 " 



5/2 U 

This asymptotics is applicable as W 3> (t — t*). The function U(t\,t) is periodic with 



respect to t\. It means the asymptotic ( pq ) is applicable on some segments of the interval 
e 4 / 5 _ <^ i. The argument of the function Wo in the neighborhood of some point 
t k as t — ► i* + is: 



A. 

e5' 



gfa) + S'fa)(t - t fc ) + O (S"(t k )(t - t k ) 2 ) 
S'(t k ) + 0(S"(t k )(t-t k )) 

~ h <bfc + C 



where 5 fe = e - 1 S'(i fc )/S'(i fc ). 

It is easy to see the argument of the function Wo may be represented as 

where is some constants depending on Sq and number k. One can see the main term 
of the asymptotics (|56| ) coincides up to shift S k with the main term of the second inner 



asymptotic expansion which is the function w(6 k ). It is easy to see for full definition of 
o 

the function U(ti,t) one must find the phase shift So. 

o 

We defined the additional constant So by matching the functions U(t±,t) and the first 

inner asymptotic expansion. 

The formula (|56|) is suitable when |Wo(#)| 5 ^ 2 » \t — When Wq{0) is small, we 

o 

consider other asymptotic formula for the function U(ti,t): 

U(ti,t) =u* + y/t=UV + <j)(t),t) . (57) 
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o 

Substitute this formula to the second-order equation for the function U{t\,t) (p8|) . 
Expand the function V (S{t)/e + 4>(t),t) with respect to the small parameter (t — t*). In 
a result the equation for the main term of the asymptotic expansion is 

p" + 6u*p 2 + = 0. 

This equation coincides with the equation for the asymptotics of the first correction of 
the first inner asymptotic expansion. The boundary conditions for the function p(S(t)/e + 
4>{t)) is obtained from the condition of the matching ( |57| ) with the asymptotics of the 
expansion (]56| ) as \9\ — ► oo. The additional constant So in the formula (|54| ) is finally 
defined at the matching of the asymptotic expansions (56), ( |5"7| ) with asymptotics of the 
inner asymptotic expansions. This get: Sq = 0. 



7 Open problems 

In this work the bifurcation of the slowly varying equilibrium of the Painleve-2 equa- 
tion was studied by matching method on the formal asymptotic approach. However it is 
necessary to note two important problems which remind out of side of our analysis. 

1. The phase shift of the oscillating asymptotic solution is undefined in our approach. 
Its definition demands much more thin calculations for the corrections of the asymp- 
totic formulas. 

2. A problem of a justification of the remainder term for the constructed asymptotic 
solution remains open as well. 
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